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Resumen. El presente documento expone algunas de las herramientas básicas desarrol¬ 
ladas para el análisis de problemas de elasticidad lineal en dos y tres dimensiones, y trans¬ 
ferencia de calor y ujo potencial bidimensionales en el entorno del programa Matlab® y 
bajo la metodología de trabajo del programa SAT-Lab® (Structural Analysis Toolbox), de¬ 
sarrollado para la enseñanza de análisis estructural y vibraciones. Es por ello que el 
énfasis del desarrollo se ha hecho en la enseñanza de los aspectos computacionales de la 
solución de problemas lineales por la vía de Elementos Finitos. Este documento consti¬ 
tuye un resumen del trabajo nal de carrera de Ingeniería de Marcos Vanella, dirigido por 
el Dr. Ing. José Inaudi, cuyo manuscrito puede obtenerse en formato pdf en la página 
www.sat-lab.com. 
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1 INTRODU CTION 

El presente documento resume el desarrollo de herramientas computacionales fe, basadas 
en el método de elementos finitos, para la solución de problemas de elasticidad lineal 
en dos y tres dimensiones, como así también para la solución de algunos problemas que 
involucran la ecuación de Laplace en dos dimensiones, particularmente transferencia de 
calor y flujo potencial. El objetivo principal perseguido en su realización ha sido facilitar 
el aprendizaje de aspectos básicos de programación del método de elementos finitos en 
la solución de problemas relativamente sencillos. Se apunta a que el usuario defina la 
secuencia de operaciones que requiere para obtener dicha solución. La utilización del 
mismo requiere del conocimiento básico del entorno y de la manera en que se introducen 
los datos y funciones disponibles, como así también de las bases de la solución de los 
problemas lineales mencionados por elementos finitos. 

El entorno del trabajo es el programa Matlab®, que básicamente es un programa 
interactivo y lenguaje de programación para realizar cálculos numéricos con vectores y 
matrices. El mismo, es actualmente, de gran utilización tanto a nivel académico como 
profesional en ambientes de ingeniería. Como las operaciones con vectores y matrices son 
partes esenciales en la solución por MEF, Matlab® resulta una plataforma conveniente 
para este desarrollo. El modo interactivo de Matlab®, reduce la velocidad de cómputo 
en algunos casos, por lo que el desarrollo no ha sido pensado para fines que requieren gran 
tarea de cálculo (mallas de decenas de miles de elementos o mayores). 

La arquitectura del programa sigue la filosofía del programa SAT-Lab®, cual es que 
el alumno lidera la solución de problemas de análisis estructural, debiendo trabajar con 
las funciones disponibles, definiendo de manera adecuada los datos de su problema, y ar¬ 
mando su camino hacia la solución. SAT-Lab® es utilizado actualmente en la enseñanza 
de análisis estructural en cursos de grado y postgrado de la Facultad de Ciencias Exac¬ 
tas Físicas y Naturales (Universidad Nacional de Córdoba) y del Instituto Universitario 
Aeronáutico. Se prevé que las rutinas desarrolladas se incorporarán a dicha plataforma 
en el futuro. 

Tanto en elasticidad lineal, como en problemas que involucran la ecuación de Laplace, 
se discretiza el dominio via elementos finitos y se transforman las ecuaciones diferenciales 
de gobierno en ecuaciones algebraicas utilizando el conocido método de Galerkin 2,3,4 de 
residuos ponderados. 

Este opera con expresiones integrales que contienen las ecuaciones diferenciales de 
gobierno del problema. Para ello se busca una solución aproximada u, que es generalmente 
un polinomio y contiene coeficientes indeterminados a ,, y que es admisible 1,2 . Para obtener 
una solución aproximada se deben determinar los valores de los a¿ de manera que el campo 
aproximado y el campo solución sean lo mas cercanos posible. Luego se procura que la 
integral de los residuos tanto en el dominio como en el contorno ponderada según una 
función de peso apropiada sea nula: 
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í Wí*R q dtt + f Wi * Rr dT = O (1) 

Jn J r 

donde W] = función de peso, Rc> = residuo en el dominio y Rr = residuo en el contorno. 
ílyT denotan dominio y contorno respectivamente. 

Por residuo se entiende la expresión de la ecuación diferencial (tanto en el dominio 
como en el contorno) valuada en la función de aproximación. Las funciones de peso 
en el método de Galerkin son las mismas que se utilizan para interpolar los valores de 
la variable incógnita, dando cada una, formulación débil 12 mediante, a una ecuación 
algebraica. Así, se obtienen integrales que corresponden a matrices y vectores globales 
del sistema algebraico, y que resultan en la práctica del ensamblaje de matrices y vectores 
correspondientes a los elementos de la discretización. 

El problema abordado por las herramientas desarrolladas es el armado del modelo 
matemático (estático o dinámico) a partir del modelo discreto que debe ser definido en 
un archivo .m de Matlab®, y el posterior análisis de dicho modelo matemático. Tanto 
modelo discreto como resultados del análisis pueden evaluarse a partir de funciones de 
graficación generadas a tal fin. 

Las rutinas fe no constituyen un programa compilado, y por ello se requiere inevitable¬ 
mente del soporte de Matlab. Incluso en casi la mayoría de las funciones generadas se 
hace uso de funciones internas de Matlab®. No obstante, es esta característica de Mat¬ 
lab, la que hace tan interesante y productivo su uso, tanto en procesos de cálculo como 
en visualización de resultados. 

Los datos relativos al modelo discreto deben ser escritos en formato adecuado al 
comienzo del archivo .m, incluyendo: Diccionario de elementos, estructura de propiedades, 
matriz de coordenadas nodales, arreglo de elementos incluyendo vector de conectividades, 
tipo de elemento, tipo de propiedad, tipo de función de cálculo de gradiente, etc., para 
cada elemento de la malla. Existen algunas funciones sencillas para importar mallas dibu¬ 
jadas en CAD, lo cual simplifica el trabajo en mallas de pocos elementos; aumentar en 
uno el grado de los polinomios de interpolación (de uno a dos), generar datos de geometría 
de elementos de contorno (sólo en problemas que involucran la ecuación de Laplace), etc. 
También dentro de lo que hace a graficación, se desarrollaron funciones que permiten 
visualizar mallas en dos y tres dimensiones, con varias opciones de gráfico. El desarrollo 
permite en los casos de estudio obtener gradientes siguiendo la lógica de programación 
anteriormente planteada y graficar el campo de soluciones y el campo de gradientes en 
gráficos de mapa de colores a través de diferentes funciones. 

2 IMPLEMENTACIÓN COMPUTACIONAL EN PROBLEMAS DE ELAS¬ 
TICIDAD LINEAL 

Es sabido que se obtiene la misma ecuación correspondiente a la forma débil del prob¬ 
lema de elasticidad lineal aplicando el método de Galerkin a las ecuaciones diferenciales 
de equilibrio dinámico o haciendo uso del principio de trabajos virtuales 1,2 para fuerzas 
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internas y externas al sistema. Dicha forma débil resulta para un elemento genérico: 


= <5{u e } 5 


[B] T {<x}.c/D e + / P(Q) [N] T [N]r/D e {ÜJ 

. ./Q e 

[N] T {F s(Q , í) }r/D e + f [N] T {f (r , í} }.dr e 


Kd(Q) [N] T [N]r/D e {ú e } 

( 2 ) 


siendo [N] la matriz de funciones de interpolación, [B] la matriz de derivadas espaciales 
de las funciones de interpolación, {er} es el vector de tensiones (tensor escrito en forma 
vectorial), es la densidad en el dominio, Kd(n) es el parámetro de amortiguamiento 
clásico en el dominio, (F 9 (q ;Í )} es el vector de fuerzas de campo distribuidas en el dominio 
y {f(r,t)} corresponde al vector de tracciones sobre el contorno T del problema. 

Aquí {u e }, {ú e } y {ü e } son los desplazamientos, velocidades y aceleraciones de los 
nodos del elemento respectivamente. En caso de haberse utilizado el Principio de Trabajos 
Virtuales para obtener (2), los <5{u e } son los desplazamientos virtuales nodales, que son 
arbitrarios, y se cancelan en la ecuación anterior obteniéndose las ecuaciones de equilibrio 
dinámico semidiscretas para el elemento. 

En el caso del presente trabajo se considera al material linealmente elástico, es decir 
{<x} = [C] {e}. donde [C] es la matriz constitutiva (tensor de cuarto orden expresado 
como matriz) del material y {s} = [B]{u e } L2 . Así, la ecuación (2) queda: 


donde, se tiene: 


mj{üj + [cj{új + [kj{uj = {rn e } + {r r J 


( 3 ) 


[™J=L,/WN] T [N]dfi« 
[cJ=4k«[N] t [N]<¡SÍ, 
M = fa e [B ] T [C] [B] 

{rne} = / n .[N] T {F í(a , t) }<i!í s 
{T«} = /r«[N] T {f (r , (| }.(¡r e 


Matriz de masa consistente del elemento. 
Matriz de amortiguamiento del elemento. 
Matriz de rigidez del elemento. 

Vector de fuerzas de campo sobre el elemento. 
Vector de cargas distribuidas sobre el contorno 
del elemento. 


El caso particular de equilibrio estático se obtiene de la ec. (3) haciendo nulos los 
vectores de velocidad y aceleración, y considerando las cargas exteriores y másicas in¬ 
dependientes del tiempo. Las matrices y vectores globales del sistema se obtienen del 
ensamble de las matrices y vectores elementales respetando la numeración global de gra¬ 
dos de libertad. 

Para el cálculo de las integrales de las matrices y vectores elementales tanto para el 
caso de triángulos, cuadriláteros (caso de elasticidad plana y axisimétrica) como hexaedros 
(caso tridimensional) se utilizan las cuadraturas Gaussianas 1 - 2,3 correspondientes. Dichos 
cálculos elementales se encuentran separados en rutinas bien definidas del programa. Por 
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ejemplo, la función de fe encargada del cálculo de la matriz de rigidez de un elemento 
cuadrilátero de 4 a 9 nodos es fequad, y está en la carpeta //e-elements del programa. 
Para el caso axisimétrico se tiene feaxquad, y en caso de triángulos de 3 a 6 nodos 
fetriang y feaxtriang; para el cálculo de matrices de masa femquad, femaxquad, 
femtriang, etc. Las matrices globales del sistema se obtienen a partir de las matrices 
elementales a través de la rutina festkcm. 

El armado del término independiente se produce de manera similar, con la diferencia 
que en el caso de cargas distribuidas sobre el contorno se tratan como actuantes sobre 
elementos del contorno que son simplemente entidades con las conectividades de los lados 
sobre el borde de elementos interiores. Una vez armado el modelo matemático, para el 
análisis del mismo pueden utilizarse funciones de SAT-Lab®. Entre los tipos de análisis se 
cuentan, solución del caso estático, obtención de modos y frecuencias naturales, respuesta 
en frecuencia, integración en el tiempo por método de Newmark, espacio de estado, etc. 
Una vez obtenidos los desplazamientos solución puede efectuarse el cálculo y alisado de 
tensiones en nodos esquina de elementos, obtención de parámetros a partir de criterios de 
falla comunes, etc. Algunos resultados gráficos obtenidos con el programa pueden verse 
en los ejemplos numéricos. Los archivos de desarrollo de dichos problemas y la manera 
de trabajo en algunos otros pueden verse en la Ref.: [1]. 

3 IMPLEMENTACIÓN COMPUTACIONAL EN PROBLEMAS DE TRANS¬ 
FERENCIA DE CALOR Y FLUJO POTENCIAL 

Resulta común que diferentes problemas físicos sean gobernados por ecuaciones diferen¬ 
ciales similares, donde varían los valores de los parámetros de las mismas, como también 
los tipos y valores de condiciones de contorno e iniciales para cada problema. Tal es el caso 
de algunos problemas de una sola variable dependiente como son transferencia de calor 
y flujo irrotacional de un fluido ideal, entre otros. Sobre dichas ecuaciones de gobierno, 
puede aplicarse el método de residuos ponderados y debilitarse las ecuaciones integrales 
mediante integración por partes 1,2,3 . La forma débil para el problema de transferencia de 
calor estacionario, considerando condición de contorno convectiva en parte del contorno 
resulta: 


V r ukVT(¿0 + / vhTdT = / vQdÜ + / vq B dV + / vhT f dT (4) 


■ ic 


- <1B 


- QC 


donde T es el campo de temperaturas a conocer, v representa una función de peso sufi¬ 
cientemente continua, [k] es la matriz de conductividades térmicas del material, h es el 
coeficiente de transmisión de calor por convección en el contorno (coeficiente de película), 
O son fuentes de calor distribuidas en el dominio, q B son flujos de calor prescriptos a 
través del contorno y T/ es la temperatura del ambiente exterior en zonas del contorno 
con convección. Suponiendo para un elemento que la temperatura en el interior del mismo 
puede interpolarse a partir de sus valores nodales según T(x,y ) = [N]{T e }, donde [N] es 
la matriz de funciones de interpolación, y aplicando el método de Galerkin la ecuación 


1507 




ENIEF 2003 - XIII Congreso sobre Métodos Numéricos y sus Aplicaciones 


(4) se transforma en el siguiente sistema de ecuaciones algebraicas: 

([ k e] + [h e ]){T e } = {r Qe }+{r gB e} + {r qhe } (5) 

donde, 

[k e ] = f 0 [B] r k[B]c/U e Matriz de conductividad térmica del elemento. 

[hj = J r [N] T /i[N]rir e Matriz de condición convectiva sobre el contorno. 

i r Qe} = jT n [N] T Qdíl e Vector de fuentes distribuidas sobre el elemento. 

{'«.} =/r..M J qBclV e Vector de condiciones de flujo prescripto. 

^Be 

{ r qhe} = f r [N] T hTfdT e Vector de condición convectiva sobre el contorno. 

Las rutinas de fe desarrolladas hasta el momento, permiten trabajar los casos tridi¬ 
mensionales (plano y axisimétrico) del problema de transferencia de calor lineal. Al igual 
que en el caso de elasticidad lineal existen funciones particulares del programa que al ser 
invocadas con los datos correspondientes calculan las matrices y vectores elementales del 
problema, integrando por la cuadratura Gaussiana adecuada. En el caso de condiciones 
de contorno dichas integrales se realizan sobre los segmentos del contorno que correspon¬ 
dan, que también son tratados como elementos con sus conectividades, tipo de elemento 
y propiedades del mismo. Las funciones ensambladoras de matrices y vectores elemen¬ 
tales (felk, felqe2nflx y felnfl2f) producen el modelo matemático resultante. Una vez 
obtenidos los resultados de la solución del caso estacionario, se calculan y alisan gradientes 
por funciones generadas a tal fin, pudiéndose luego graficar dichos resultados. El trabajo 
en casos de transferencia de calor puede verse en ejemplos de la Ref.: [1]. Se muestran 
resultados en un caso concreto en la sección de ejemplos numéricos. 

En el caso de flujo potencial estacionario, considerando función potencial de veloci¬ 
dades, la forma débil de su ecuación de gobierno resulta: 

f V T vI V(¡)dü= í vV 0 dr (6) 

donde v es la función de peso, [I] es una matriz identidad, ó es la función potencial de 
velocidades (campo escalar) y V 0 es la velocidad del flujo normal a la parte del contorno con 
condición de contorno secundaria. La aplicación del método de Galerkin a la anterior, 
considerando que (/>(x,y) = [N]{0 e }, lleva al siguiente sistema de ecuaciones para el 
elemento: 

[k e ]{0 e } = {r qBe } (7) 

donde se tienen: 

[k e ] = f Q [B] r I[B]c/U e Matriz de difusividad del elemento. 

{ r q ge } = f r [N] r Vodr e Vector de condición de velocidad normal prescripta. 
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Un razonamiento análogo se puede efectuar para el caso de trabajarse el problema 
de flujo potencial plano mediante la función de corriente (campo escalar ortogonal a la 
función potencial). En el caso axisimétrico sólo se ha desarrollado el trabajo con función 
potencial ya que la función de corriente no satisface la ecuación de Laplace en coordenadas 
cilindricas. Las rutinas de armado de matriz y vector elementales ([k e ] y {r qBe }) resultan 
las mismas que en el caso de transferencia de calor, como así también las rutinas de 
solución del problema estacionario y cálculo de gradientes (en este caso las velocidades). 

Tanto en transferencia de calor como en flujo potencial, resulta conveniente a veces 
modelar dominios infinitos o seminfinitos. Tal es el caso de modelarse la transferencia de 
calor desde o hacia un objeto enterrado en el suelo, flujo alrededor de un perfil en vuelo 
libre, etc. Dicha condición de dominio infnit-o en ciertas direcciones puede modelarse 
a través de elementos infinitos 1,3 . En el programa fe se han desarrollado las rutinas 
correspondientes para poder efectuar el análisis estacionario sobre este tipo de modelos. 
La implementación detallada en este tipo de casos puede verse en la Ref.: [1]. Los 
resultados obtenidos para un ejemplo particular pueden verse en la sección de ejemplos 
numéricos. 

4 CONCLUSIONES 

La meta propuesta como objetivo del presente trabajo, de generar un paquete de her¬ 
ramientas computacionales de trabajo y aprendizaje para los problemas estudiados ha 
sido alcanzada. La lógica de programación desarrollada probablemente diste mucho de ser 
la óptima en cuanto a "tiempo de máquina", pero es sencilla y permite su entendimiento 
a quien desea conocer los procesos básicos de cálculo de subrutinas que indefectiblemente 
se encuentran en un programa de Elementos Finitos, lo cual fue uno de los objetivos 
fundamentales de este trabajo. 

Quedan muchas aristas por donde avanzar en el desarrollo, la intención es que el avance 
y la tarea programación sean continuos, para poder anexar a SAT-Lab® en el corto plazo 
las carpetas con herramientas relativas a elementos finitos para problemas de elasticidad 
lineal y no lineal. 

5 EJEMPLOS NUMERICOS 

5.1 Ejemplo 1: Estudio de la concentración de Tensiones sobre una llave de 
ajuste cargada. 

El presente ejemplo es un ejemplo típico de elasticidad lineal tridimensional en tensión 
plana. Consiste en una llave de ajuste mostrada en figura (1), la cual se encuentra 
sometida a carga dirtribuida normal q n = 0.3 [kgf/mm 2 ] sobre una longitud de 30 [mm] 
en el contorno superior del mango, tal como se muestra en dicha figura. El espesor de la 
llave es t = 10 [mm]. El material con que está construida la llave es Acero SAE 1040. Se 
pretende analizar el estado tensional de la llave bajo la carga mencionada. 

En el modelo se supondrá que los desplazamientos en las dos caras en contacto con 
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Figura 1: Ejemplo 1: Llave de ajuste bajo carga. 

la tuerca están restringidos a un valor nulo, suposición que no es del todo cierta ya 
que la distancia entre dichas caras se mantendrá constante y por ello no se registrará 
concentración de tensiones en el punto central de dicho contorno parcialmente hexagonal. 
Para considerar condiciones de contorno principales más realistas deberían tenerse en 
cuenta restricciones nodales multipunto. 

El archivo de datos de este ejemplo puede verse en el Apéndice [A] de la Ref.: [1]. 
Se modela el dominio con elementos triángulos de seis nodos. Se definen diccionario de 
elementos y propiedades, y se importa la malla de elementos de un archivo de CAD, la 
cual puede verse en la figura (2). 


604 nodes, 259 elements 



Figura 2: Ejemplo 1: Malla de elementos triangulares. 

Luego se define la matriz RESTRAINTS con las restricciones en grados de libertad 
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de nodos de los contornos en contacto con la tuerca, y se genera la matriz DOFS de 
numeración global de grados de libertad. Posteriormente se calcula la matriz de rigidez 
del sistema. 

Se obtienen los elementos que pertenecen al contorno mediante la función febelenrents- 
gen, y se determina sobre cuales segmentos del contorno está aplicada la carga. Así se 
calcula el término independiente con las cargas nodales equivalentes actuantes sobre dichos 
elementos del contorno. 

Luego se soluciona el sistema de ecuaciones y se calculan y alisan tensiones. Por último 
se gráfica en este caso la tensión de comparación obtenida por criterio de Huber-Von 
Mises respecto de la tensión de fluencia del material (inversa del coeficiente de seguridad). 
Pueden verse claramente los puntos donde hay concentración de tensiones. Dicho gráfico 
se presenta en la figura (3). 



Figura 3: Ejemplo 1: Tensión de Comparación según criterio de Von Mises respecto de tensión de fluencia. 


Para disminuir las concentraciones de tensión en los puntos críticos a ambos lados del 
mango deberían suavizarse los quiebres que sufre el contorno en dichos puntos, que a 
pesar de ser muy leves, tienen una gran influencia sobre el comportamiento tensional en 
la sección que los une. El gráfico de la estructura deformada, cuyos desplazamientos se 
han amplificado 50 veces puede verse en la figura (4). 
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604 nodes, 259 elements 



Figura 4: Ejemplo 1: Estructura deformada (desplazamientos amplificados 50 veces). 

5.2 Ejemplo 2: Estudio de la pérdida de calor de tubo enterrado. 

El presente ejemplo ha sido tomado de Ref.: [9] donde se ha resuelto utilizando factores de 
forma para la conducción. Se tiene un tubo de 6 pulgadas de diámetro exterior enterrado 
con su eje central 30 pulgadas abajo de la superficie del suelo que tiene una conductividad 
térmica promedio k = 0.20 [Btu¡(hr.pie.°F)\. Se pretende determinar la pérdida de calor 
por pie de longitud de tubo, si la temperatura de la superficie del tubo es de 200 °F y la 
temperatura de la superficie del suelo es de 40 °F. 

El proceso de solución, utilizando fe, requiere de la utilización de elementos infinitos ya 
que el suelo en profundidad y extensión se modela como un dominio infinito. El problema 
resulta simétrico respecto de un plano que pasa por el eje del tubo y es perpendicular 
a superficie del suelo. Luego se modela una mitad del dominio. El dominio modelado 
utiliza elementos triángulos y cuadriláteros, y además elementos infinitos. El código de la 
solución de este ejemplo se expone en el Apéndice [B] de la Ref.: [1]. 

Resumidamente, se comienza definiendo propiedades de elementos interiores y se im¬ 
porta la malla inicial. Se obtienen los elementos del contorno iniciales y se genera el vector 
de grados de libertad inicial. Luego se define cuales de los elementos del contorno serán 
lados de elementos infinitos, se define como polo u origen de decaimiento 1,3 de algunos 
elementos infinitos al punto situado en el eje del tubo, y de otros, puntos tales que el lado 
del útlimo elemento infinito sea la línea del horizonte del suelo. Luego se actualizan los 
datos teniendo en cuenta dichos elementos. La malla resultante puede verse en figura (5). 

Posteriormente se calcula la matriz de coeficientes del sistema, se restringen los nodos 
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Figura 5: Ejemplo 2: Malla teniendo en cuenta elementos infinitos. 

de lados con temperaturas prescriptas y se soluciona el sistema de ecuaciones a través de 
la función sesolve de Sat-Lab®. Vemos que en este caso el término independiente tiene 
todos sus valores nulos ya que no hay condiciones de flujo no nulas en el contorno ni 
fuentes internas. 

Obtenido el campo de temperaturas resultante se calculan y alisan gradientes. Se 
integran numéricamente los gradientes en nodos esquina de elementos que pertenecen al 
contorno del tubo y duplicando este resultado se obtiene la pérdida de calor por pie de 
longitud de tubo: 

>>Matlab 

Flujo obtenido en ejemplo de libro Kreith= 

66.5 BTU/(hr*pie) 

Flujo obtenido por elementos finitos= 

66.683 BTU/(hr*pie) 

Vemos que el valor obtenido concuerda con el resultado dado en Ref. [9]. El campo de 
temperaturas para este ejemplo puede verse en la figura (6). 
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Figura 6: Ejemplo 2: Campo de temperaturas resultante. 

5.3 Ejemplo 3: Caso axisimétrico: Modelo de flujo potencial para Túnel de 
viento de Baja Turbulencia. 

En el presente ejemplo se desea tener un conocimiento aproximado de las velocidades 
en cámara de ensayos de un túnel de viento de baja turbulencia de cámara de ensayos 
cerrada y relación de contracción en el convergente de 36 a 1. El plano corresponde 
al túnel de baja turbulencia del Centro de Investigaciones Aplicadas, dependiente del 
Instituto Universitario Aeronáutico. La descripción de túneles de viento de este tipo 
puede verse en Reís.: [10] y [11]. 

En el modelo simplificado del siguente ejemplo se considera flujo potencial, además 
que no existen pérdidas de presión en mallas de la cámara de tranquilización y que la 
cámara de ensayos se encuentra vacía. El modelo resulta simétrico respecto al eje del 
túnel y se tratará con elementos triángulos axi s i m étricos utilizando la función potencial 
de velocidades. 

Las condiciones de contorno serán las que se ven en la figura (7). Se tendrá flujo normal 
nulo sobre el eje se simetría y sobre la línea que representa las paredes del túnel. Además 
se restringe la función potencial a un valor nulo en los nodos sobre la salida del divergente 
y se especifica una velocidad normal de Uo = 1250 [mm/seg\ en la sección de entrada al 
túnel. 
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Figura 7: Ejemplo 3: Modelo de Túnel de Baja Turbulencia y condiciones de contorno. 

El trabajo con fe sigue la lógica vista, con la excepción que se definen los elementos 
como axisimétricos, pudiéndose ver el código de la solución en el Apéndice [B] de la 
Ref.: [1]. La malla importada de elementos triangulares puede verse en la figura (8) y 
las velocidades calculadas en la figura (9). La velocidad máxima obtenida en cámara de 
ensayos resulta de módulo: V max = 47.26 [m/seg]. 


3636 nodes, 1683 elements 



Figura 8: Ejemplo 3: Túnel de Baja Turbulencia. Malla de elementos triangulares axisimétricos (medidas 
en mm). 


Además, suponiendo que el fluido tiene densidad del aire p = 1.225 [ kg/m 3 ] (atmósfera 
standard), puede graficarse la caída de presión estática a lo largo del túnel, que es de 
soplante aspirante conectada a la salida del difusor, siendo esta caída de presión igual a 
la presión dinámica que se tiene en cada punto (modelo de flujo potencial). Dicho gráfico 
puede verse en la figura (10) resultando la caída de presión en el punto de cámara de 
ensayos con mayor velocidad del fluido A P est = 1368.02 [ Pa ]. 

La comparación de los resultados obtenidos con los resultados de un modelo de flujo tur¬ 
bulento y condición de deslizamiento no nulo en las paredes, desarrollado en el programa 
CFDesign®, muestra que, aunque las hipótesis de flujo son diferentes, los resultados son 
similares en este caso. Dicho modelo arroja una velocidad máxima en cámara de ensayos 
de 48.5 [m/seg] aproximadamente, tal como se muestra en el gráfico de CFD Desing® 
en la Ref.: [1]. Este valor resulta levemente superior al obtenido por flujo potencial. 
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Figura 9: Ejemplo 3: Túnel de Baja Turbulencia. Módulo de velocidades obtenidas. 
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Figura 10: Ejemplo 3: Túnel de Baja Turbulencia. Caída de presiones estáticas. 
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